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In this paper we study the energy spectrum of a two dimensional electron gas (2DEG) 
in a two dimensional periodic magnetic field. Both a square magnetic lattice and a trian- 



^ . gular one are considered. We consider the general case where the magnetic field in a cell 

cn ■ 

OO ' can be of any shape. A general feature of the band structure is bandwidth oscillation as 

O' 

t^^ ' a function of the Landau index. A triangular magnetic lattice on a 2DEG can be realized 
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^ . by the vortex lattice of a superconductor film coated on top of a heterojunction. Our 



calculation indicates a way of relating the energy spectrum of the 2DEG to the vortex 
Vh ' structure. We have also derived conditions under which the effects of a weak magnetic 

^ . modulation, periodic or not, may be reproduced by an electric potential modulation, and 
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Q , vice versa 

> 

X 

^ : PACS: 71.25-s, 73.20DX, 74.60.-w 

To be published in Phys. Rev. B. 



I. INTRODUCTION 

The behavior of a 2DEG in a homogeneous magnetic field modulated by a periodic 
electric potential or magnetic field is a subject of some recent investigations [1-7]. In the 
case of a one dimensional electric modulation (IDEM), resonance between cyclotron radius 
Re and the period of modulation a results in a novel magnetoresistance oscillation[l,2]. It 
was found that the electrons do not move along the trench formed by the modulation po- 
tential when 2Rc = {k — l/4)a, and move with maximum velocity when 2Rc = {k + l/4)a, 
where k is an integer. Because Re changes with the magnetic field, the magnetoresis- 
tance oscillates when the strength of the magnetic field increases. This result has been 
explained quantum mechanically by Gerhardts et a/[l], and classically (when Re ^ a) 
by Beenakker[2]. A variant system that uses a one dimensional magnetic modulation 
(IDMM) instead of a IDEM shows similar behaviors[3]. In this case the electron drift is 
absent when 2Rc = {k + l/4)a, and a maximum drift velocity along the trench is obtained 
when 2Rc = {k - l/4)a. 

A delicate structure in the electron energy spectrum emerges when the modulation is 
two dimensional. It is well-known that the spectrum in a magnetic field with a two di- 
mensional electric modulation (2DEM) exhibits fractal behavior — the so-called Hofstadter 
spectrum[4]. The fragmentation of bands will suppress the part of magnetoresistance that 
comes from the variation of band conductivity, and only a weaker, diffusion-related mag- 
netoresistance oscillation can be observed[5]. Fractal spectral structure also exists for a 
two dimensional magnetic modulation (2DMM), which was studied by Wu and UUoa in a 
recent study. They also studied the collective excitations of the electron gas, and found 
that it does not map out the Hofstadter spectrum exactly [7]. Although it would be inter- 
esting to observe the excitation spectrum of a 2DMM system, there is still no attempt on 
an experimental realization up to the writing of this paper. 

In this paper we study the energy spectrum of the 2DEG as a function of quasi- 
momentum, and show that there is a one to one correspondence between it and the spatial 
distribution of the magnetic field. Both a rectangular magnetic lattice and a triangular 
one are considered. The latter can be generated using a superconductor in a vortex state. 



Therefore we can have a periodic modulation on a 2DEG by coating a superconductor 
film on top of a semiconductor heterojunction. Using the connection between the energy 
spectrum and the shape of the magnetic modulation, we can use the 2DEG as a probe for 
getting some basic parameters of the superconductor film. 

We also explore the possible connection between electric and magnetic modulations. 
There is an exact mapping at the operator level between the Hamiltonians of a IDMM 
system and a IDEM system[9]. No such simple mapping exists between 2D modulation 
systems. However, the matrix elements of the Hamiltonian of a 2DMM are similar in 
structure to those of a 2DEM. Hence, it is possible to have a formal connection between 
these two, albeit in a nontrivial way. The cases with one dimensional, two dimensional 
periodic modulations, and non-periodic field distributions are discussed. 

This paper is organized as follows: In Sec. 2 we study systems with weak perturbations. 
In Sec. 3 the exact energy spectrum is obtained by diagonalizing the Hamiltonian matrix. 
In Sec. 4 a way of relating the vortex structure to the behavior of a 2DEG is discussed. In 
Sec. 5 we discuss the correspondence between the electric modulations and magnetic field 
modulations. The last section is a conclusion. 

II. WEAK MODULATION 

A. Rectangular lattice 

Consider a 2DEG in a homogeneous magnetic field Bq perturbed by a weak rectangular 
magnetic lattice B{x,y). In order to simplify the notation, we choose magnetic length 
A = (h/eBo)^/'^ to be the unit of length, Bq to be the unit of magnetic field, and huJc to 
be the unit of energy, where Uc = {eBo/mY''^ is the cyclotron frequency. This is the most 
natural choice of units in the sense that all physical quantities are around the order of 
unity in the regime of interest. The Hamiltonian can be written as 

^=2 V 'a^'^ + ^^^^^'^V ^2 V 'a^ + ^^^'^'^V ■ ^^^ 

We have choosed the vector potential of the homogeneous magnetic field to be Ao{x, y) = 
—yi. A{x, y) = Ax{x, y)i + Ay{x, y)j is the vector potential of the perturbing field. 



An eigenfunction of the unperturbed Hamiltonian for the n-th level is 

y L 

where A^^ is a normalization constant, L is the size of the system, and ki is a constant. 
Since the perturbed system has rectangular symmetry, we would like to construct from (2) 
an unperturbed eigenstate with the same symmetry. Assume ai and a2 to be the lattice 
constants. Discrete translation operators that commute with the Hamiltonian (1) are the 
magnetic translation operators Ti, T2. 

_ . (3) 

T2 = e^"2^T2, 

where Ti and T2 are the usual translation operators. An unperturbed eigenstate that 
satisfies 

is found to be 

CXD 

^n.,.. = Y. e-^^"^"^^^,,,_-. (5) 

1^ — 00 

Eq.(4) defines the quasi-momenta ki and /c2, which are good quantum numbers. We 
assumed there is only one fiux quantum $0 = h/e per unit cell in getting Eq.(5). That 
is, we used aia2 = 2/t in the new unit system. Then, a state in a Landau level is only 
coupled to states in the other Landau levels with the same (/ci, A;2), relieving the necessity 
of degenerate perturbation. The case with arbitrary magnetic fiux per unit cell is more 
complicated, and we will comment on it later. 

The first order energy perturbation can be found by calculating the diagonal matrix 
elements of the Hamiltonian on the unperturbed basis: 
^ 00 

Lm= — 00 



Zlm = 27r 



(l)2 + (!!^)2 

tti a2 



(6) 



where B{1, m) are the Fourier components of B{x, y), 

B{x,y) = ^S(/,m)e2^*^^/"ie2^'""^/«^ (7) 

Im 

and L^ is the associated Laguerre polynomial {L^i = 0). Notice that B{x,y) is in units 
of So- Eq.(6) can be readily reduced to the result of a system with a one dimensional 
sinusoidal perturbation in Ref. 3. 

Using the asymptotic form of the associated Laguerre polynomials in the large n limit, 
it can be shown that, in the 4-point approximation (where all the Fourier components are 
zero except -B(±l, 0) = -B(0, ±1) 7^ 0), the bandwidth is zero when l^fmi = (k + l/4)7r, 
and reaches maximum when I^Jrm = (k — l/4)n. Therefore the widths of the energy 
bands in Eq.(6) oscillate with the Landau index. Physically, it is clearer to write these two 
conditions as 2Rc = {k± l/4)a, where the radius of cyclotron motion R^ — v2n when n 
is large, and a = v27r. These are the same as the resonant conditions for a IDMM. The 
energy bands will be further split into p sub-bands when there are a rational number (say 
p/q) of flux quanta in a unit cell. A proper unperturbed wave function for such a system 

is 

p 00 

j = l l= — oo 

The coefficients dj are to be determined from a Harper-like equation after we turn on the 
perturbation[10], more details can be found in Ref. 7. 

B. Triangular lattice 

An easy way to realize a two dimensional periodic magnetic modulation would be 
using a superconductor vortex array with the symmetry of a triangular lattice. We choose 
Ao(x, y) = —yt to make our Hamiltonian explicitly translational invariant in x. The 
magnetic translation operators that commute with the Hamiltonian of this system are 

Ti = giTi 



T2 = e^^""/'(72T2, 



(9) 



where gi = exp(— ztt), and g2 = exp(— i7r/2) are gauge factors for shifting the center of a 
hexagonal unit in Sn{ki, /C2) to the origin. This will become clearer in a later discussion. 
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Instead of putting gi and §2 in the magnetic translation operators, we can put them on 
the right hand side of Eq.(4). This amounts to a redefinition of the quasi-momenta. 
An unperturbed eigenstate with triangular symmetry is given by 

oo 

^nk,k, = Y. (-l)^^^-'^/'e-('=^/2-^^)>„,,^_a^, (10) 

again we assumed one flux quantum per unit cell (v3a^/2 = 27r). The case of one-half 
flux quantum per unit cell will be discussed at the end of this section. 
The flrst order energy perturbation is 

^ oo 

r,s^ — oo 

/ 2 2\ 

Wrs = —]=\^ - sr + r ) 
v3 

B(r, s) are the Fourier components of B{x^y), 

B{x,y) = J2Bir, s)e*(->2/i)-('^f''i+-&2)^ (12) 

r,s 

where bi = {2ii /a){i — j/\/?>) and 62 = (27r/a)(2j/v^) are unit vectors of the reciprocal 
lattice. 

If we let all of the Fourier components be zero except six of them (we will call this 
a 6-point approximation): -B(±l, 0), -B(0, ±1), and i?(±l, ±1), which are all equivalent to 
each other. Then we have 



27r 27r 27r 

2cos( — x) cos( ^ y) + cos( _ 2-0 



(13) 



i?(x,2/) = 2S(l,0) 
The corresponding energy perturbation is 

Sn{ki,k2) = -e-"/^S(l, 0){Ll + Ll_^)[cos{aki) + cos{ak2) + cos{a{ki - ^2))]. (14) 

Notice that ki and /c2 are along non-orthogonal directions Ti and T2, we need to make a 
transformation 



X 



ki — k. 

, kx v^, (15) 

k2 = Y + -Y^y 



to put £n on an orthogonal basis k^ and ky. See Fig.l for plots of B{x, y) and £^0(^x5 ky), 
where -B(l, 0) = 0.05. Notice that the magnitude of variation of B{x, y) is 95(1, 0) (from 
-3S(1, 0) to 6S(1, 0)), while that of ^o(A:x, ky) is 4.5e-^/^S(l, 0), which is a much smaller 
variation. It appears that the electrons sense the average magnetic field more than they 
sense the fiuctuation. The plot in £o{kx, ky) will be shifted to the right by n if both factors 
gi and (72 in Eq.(9) are identity. 

In a real physical situation, the triangular array is formed by a superconductor vortex 
state, where the magnetic fiux per plaquette is $o/2 instead of $o- In this case the 
magnetic translation operators defined in Eq.(3) do not commute (T1T2 = — T2T1), and 
the quasi- momenta in Eq.(4) are no longer good quantum numbers. We have to choose 
a unit cell that consists of two magnetic plaquettes to have a mutually commuting set of 
H, Ti and T2. Assume a is the lattice constant of the magnetic plaquette, then there is one 
flux quantum $0 within the area {V3/2)aia2 of the unit cell, where ai = 2a, 02 = a. The 
energy perturbation for this system is similar to Eq.(ll), but with two differences: Firstly, 
(—1)^* is replaced by (—1)^^'', and thus can be dropped. Secondly, kia in the exponent 
is replaced by /ciai, and /c2« by 2/^202- Hence there are two identical regions in the first 
Brillouin zone. This is related to the following fact: Instead of choosing ai = 2a, we can 
choose a2 = 2a, and ai = a. The overlapping region of these two different Brillouin zones 
have the same energy spectrum because we are free to enlarge our unit cell in either ways. 
The dispersion surfaces in the other half regions of either Brillouin zones just replicate 
the energy spectrum in the overlapping region. The existence of two identical regions in a 
Brillouin zone is related to the two-fold degeneracy of such a system. 

III. EXACT ENERGY SPECTRUM 

The exact energy spectrum for an electron in a magnetic field with strong modula- 
tion can be found by diagonalizing the Hamiltonian (1) numerically. For a square mag- 
netic lattice, the result of diagonalization in the 4-point approximation is shown in Fig. 2. 
The energies from the perturbation calculation are lower than the exact energies since a 
positive-definite contribution from the quadratic terms of the vector potential is neglected 



in calculating Eq.(6). Notice that the magnetic modulation is of the same strength as the 
homogeneous component, but the energy shift is just a fraction of huc- 

Matrix elements of the Hamiltonian (1) with triangular modulation are listed below 
for reference. Write the Hamiltonian as if = Hq + Hi + H2, where Hi is linear in A, and 
H2 is quadratic in A. The matrix elements of Hi on the unperturbed basis (10) are found 
to be 



{n'\Hi{ki,k2)\n)^^,y^^ 



\Wr 



n' + n 



Wr 



\Wrs) 



2n' 



Wr 



K-ri^rs) 



Gn'n{r,s). 



(16) 



The matrix elements of H2 are 

(n'|i72(/ci,A;2)|n) („,>„) 
= Y, [Mr, s)A,{r', s') + Ay{r, s)Ay{r', s')]L;^'""(w,+,/,,+,/)G„.^(r + r' , s + s'), C^^) 



where the function Gn'r, is defined to be 



n! 



Gn'nir.s) = \ —{V2n 



n'\ 



^r,-n, ^r ^ 2s -r 
a y/3a 



-Wrs f^ ( _-l'\rs is{kia+Tr) —ir{k2a+TT) 



(18) 



The vector potential A(r, s) is related to B{r,s) by B{r,s) = iq x A(r, s), where q 
rbi + sb2 7^ 0. Using the Coulomb gauge condition, we have 



A(r, s) 



kT 



X B{r, s). 



(19) 



Under the 6-point approximation, there is only one independent component, say Aj:{l,0), 
of the vector potential, which is related to S(1,0) by Aj;{l,0) = {\^a/8ni)B{l,0). The 
energy spectrum for a system with a 6-point approximation is shown in Fig. 3, the mod- 
ulation strength 5(1,0) is chosen such that the minimum of the total magnetic field 
B(x, y) = Bq + B{x, y) is zero. Numerical calculation shows that the transition ampli- 
tudes between different Landau levels oscillate with respect to the level separation An 
when An is small, and drop quickly to zero beyond An ~ 10. The asymptotic behav- 
ior goes roughly like l/(An)!. We were careful to choose the rank of the matrix being 
diagonalized large enough to ensure the convergence of our result. 



When the variation of magnetic field is large, adjacent energy levels begin to mix with 
each other. This can be seen from Fig. 4. The modulation part is nine times larger than 
the homogeneous part, therefore in some region of the unit cell the total magnetic field 
reverses the direction. Notice that these levels in general do not touch. This conforms to 
a theorem by Von Neumann and Wigner, which states that, it needs to vary at least three 
parameters for a generic Hamiltonian to have an accidental degeneracy [11]. 

Aharonov and Casher once showed that the ground state energy of a spin- 1/2 particle 
in a nonuniform magnetic field is zero [12] . This result is checked numerically in this periodic 
case by adding a spin interaction term {az/2)B{x,y) to the Hamiltonian (1). Since (10) 
is not a very good basis for this new Hamiltonian, it is necessary to diagonalize a larger 
Hamiltonian matrix to get a correct result. The lowest dispersion curve we get for this 
Hamiltonian with spin term is found to be fiat and zero with negligible deviation, which 
agrees very well with Aharonov and Casher's result. 

IV. ENERGY SPECTRUM AND VORTEX STRUCTURE 

It is not hard to see that, when the energy perturbation is small, we can determine 
B{x,y) with the help of Eq.(ll) and Eq.(12) if the energy spectrum Sn{ki,k2) over the 
entire first Brillouin zone can be mapped out precisely. Therefore, in principle, the 2DEG 
beneath the superconductor film can be a probe for measuring the vortex structure. How- 
ever, this is a difficult experimental task compared to other methods like neutron diffrac- 
tion, scanning tunneling microscopy[13] and electron wave holography [14]. Angle Resolved 
Photoemission Spectroscopy (ARPES) can be a tool of measuring the band structure, but 
it can barely probe beyond lOA below the surface, while our 2DEG is buried down in the 
heterojunction. Besides, the electron density of a conventional 2DEG (~ 10^^ cm~^) is 
too low for ARPES to have the precision we need. In spite of these difficulties, at least the 
bandwidth can be measured with confidence, and this can give us some information about 
the vortex state. In this section we will take a simple model of vortex structure as our 
starting point, and demonstrate how the bandwidth of the 2DEG may help us determine 
the coherence length of the superconductor. 
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A convenient choice for the structure of a vortex was proposed by Cleni[15]. The 
magnetic field distribution of a single vortex in the Clem model is 

where Kq and Ki are the modified Bessel functions, p is a radial coordinate, A is the 
penetration depth, and ^„ is a variational parameter that is of the same order as the 
coherence length. ^„ and A are the only two parameters in this model. The magnetic field 
of a vortex array is Bs{r) = '^]^b{r — R), where the summation is over all sites of the 
vortices. Therefore, the Fourier components, 

where Aceii is the area of a unit cell, q = rbi + sh2, and Q = {\q\^ + A~^)^/^. 

The magnitude of oscillation of the magnetic field Bs{x,y) at the surface of the su- 
perconductor is attenuated through the distance d between the superconductor film and 
the 2DEG. Fourier component of wave vector q is damped out by a factor exp(— |g|(i)[l]. 
Therefore, the Fourier components of the magnetic field that appear in Eq.(ll) are ac- 
tually B{r,s) = Bs{r,s) e~'^'^, where \q\ = 47r/(\/3a) Vr^ — rs + s^. For example, if 
a = 1500 A (for Bq = 0.1 Tesla), and d = 400 A, then S(l, 0) = 0.165^(1, 0), S(l, -1) = 
0.04i?s(l, —1). The modulation on the 2DEG is small unless Bq is much less than 0.1 T, 
or d is much less than 400 A. Therefore in most of the cases the perturbation formula 
works well. 

Assuming that the vortex array is dense enough so that the 6-point approximation is 
valid, then the width of the n-th band is 

ASr. = 4.5e-"/^S(l, 0){Li{2n/V3) + Li_^{2n / Vs)) . (22) 

One way of measuring the bandwidth is by measuring the differential Hall conductivity of 
the 2DEG, which in the coUisionless limit is 

^HiE) = -^D{E)f{E), (23) 
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where D{E) = ^^ / j:j^d{E — Sn{k)) is the density of states, and f{E) is the Fermi distri- 
bution. Notice that we have put back the real units in the above formula. This expression 
can be derived from the Kubo formula, under the assumption that the modulation is small 
compared to the average field. The effect of modulation appears through the S{k)'s in the 
density of states and in the Fermi distribution. Eq.(23) is the same as the classical formula 
o'H = Pee/Bo, where pe is the density of mobile electrons. 

In the presence of disorder, D{E) has to be replaced by the density of extended states 
De{E), or more precisely, the density of Chern numbers [16]. For a homogeneous magnetic 
field, it has been shown that De{E) is a delta function for each Landau levels. When the 
magnetic field is modulated, we expect that the width of De{E) be approximately equal 
to the bandwidth in the weak disorder limit. A tricky point in measuring the bandwidth 
is that we need a large B to have a good resolution in a quantum Hall system, but 
that tends to pack the vortex array tightly and shrink the bandwidths. Therefore it is 
necessary to find a compromise between these two constraints. A theoretical calculation 
of the Hall conductivity according to (23) at zero temperature is shown in Fig. 5. It varies 
continuously when the Fermi energy is inside the band, and takes quantized values when 
the Fermi energy is in the gaps. The slopes are discontinuous at the band edges since the 
density of states are discontinuous for a 2D system. How a finite disorder will change the 
extended density of states remains an open question. 

From the bandwidths we can get one Fourier component 5(1, 0) and hence an alge- 
braic relation between ^„ and A through Eq.(21) and Eq.(22). This is useful because the 
coherence length of a type-II superconductor with large k is usually more difficult to be 
measured than the penetration depth. Therefore, if A is known beforehand, measurement 
of the bandwidth can give us the coherence length. The measured values of 5(1, 0) from 
different bands can be used as a consistency check of this approach. In principle it is 
also possible to determine the temperature dependence of ^„ since the bandwidth will be 
different when ^„ is changed by temperature. 

The Hall conductivity gives more information about the band structure than just 
the bandwidths. One may compare the density of states measured according to Eq.(23) 
and that calculated from the theoretical model. If the sample is clean enough, such a 
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comparison should yield information about the higher Fourier components of the magnetic 
field, and therefore more information about the magnetic flux lattice. 

When the strength of the modulation is large, neighboring bands are mixed due to 
strong oscillations (see Fig. 4), also Eq.(23) for the Hall conductivity is no longer valid. 
However, each filled band still contributes an integer number of e^ /h to the Hall con- 
ductivity as long as it is not in touch with the other bands (usually adjacent bands are 
separated by avoided crossings). These avoided crossings (or, local gaps) can be closed by 
tuning S to a particular value B* . It does not violate the Von Neumann-Wigner theorem 
since three parameters {k^^ky^B) are used to induce an accidental degeneracy. When a lo- 
cal gap collapses, the Hall plateau that corresponds to that gap disappears. A new plateau 
emerges when the gap is re-opened by increasing B away from B* . This plateau may or 
may not stick to its original value. In case that it does not, (Jh will jump up or down 
by an integer multiple of e^//i[17]. This dramatic effect can also be achieved by using an 
electric modulation to control the band crossings. 

V. CORRESPONDENCE BETWEEN GENERAL MAGNETIC AND ELEC- 
TRIC MODULATIONS 

A. One dimensional modulation 

It is well-known that the equation of motion of a 2D electron in a constant magnetic 
field can be reduced to that of a ID electron in a parabolic electric potential. This analogy 
can be generalized when the magnetic field is modulated along one direction[9], where the 
corresponding electric potential V{x) is determined by 



ax 

For a weak IDMM with B{x) = 1 + i?cos(27r/aa:), the corresponding electric potential 

■r./ N (x — c)^ ~, . a . 27r , ^, 

V{x) ~ ^ — —^ + B{x - c) — sm — X. (25) 

2 271 a 

Besides the "all-magnetic" and "all-electric" systems, it is also possible to transform the 
term (x — c)^/2 in Eq.(25) back to a homogeneous magnetic field, and keep the second term 
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intact. This is equivalent to a mixed system with a magnetic field plus a modulating electric 
field. Following Beenakker's analysis on the drift of the guiding center of an electron in a 
1DEM[2], the mean square drift velocity in the classical limit Re ^ a is found to be 

i'^d) = (B — ) —COS ( -), (26) 

2n a a 4 

which is the same as the one in Xue's paper [3]. It has maxima under the resonant condition 
2Rc = (K-l/4)a, 

B. Two dimensional periodic modulation 

There is no simple mapping as in Eq.(24) when the modulation is bidirectional. How- 
ever, a connection still exists when the perturbation is small. Consider a 2DEG in a rect- 
angular 2DEM. The first order energy perturbation is similar to that in Eq.(6). The only 
difference is that B{l,m) [Ll^{zim) + L\_i{zim)) is replaced by 2Vi{l,m)Ln{zim)- There- 
fore, as far as the energy spectrum goes, a magnetic modulation can be simulated using 
an electric potential whose Fourier components satisfy 

2Vi{l,m)Ln{zim)=B{l,m) [Ll{zim) + Ll_^{zi„,)] (27) 

It can be readily reduced to a one dimensional relation if all the components with / ^ 
are zero. Unlike the effective potential in Eq.(25), this one is valid only in the one-band 
approximation, and is different for different bands. 

A curious feature of this formula is that if a zim happens to coincide with a zero of 
the Laguerre polynomial, the (/,m)th Fourier component of the effective potential does 
not exist, simply because there is no way a Vi{ljm) can contribute to the broadening of 
the n-th level, no matter how large it is. Therefore, for a given Landau level with n 7^ 0, 
an electric modulation can only mimic a magnetic modulation whose Fourier components 
B{l,m) vanish when zim are at zeros of L^. Conversely, a magnetic modulation can only 
mimic an electric modulation whose Fourier components Vi{l,m) vanish when zim are at 
zeros of L^ -|- i^^_i. 

C. General two dimensional field distribution 
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It is possible to extend the above connection to a random magnetic field. The only 
restriction is that the random component be weak compared to the homogeneous com- 
ponent. Consider the matrix elements of the Hamiltonian in the unperturbed basis (2) 
(no requirement on one flux quantum per unit cell) . Express the magnetic field in Fourier 
integral instead of Fourier series since there is no discrete translational symmetry here, 
then we have 

(n, q - q'\Hf{x, y)\n, d) = \ j dk2B{q\ k2){Ll + L^_i)e-^('^"+'=')e^'^''=^e^'=^^ (28) 

where the argument in the associated Laguerre polynomial is [q'"^ + /c|)/2. For a general 
electric potential V{x,y), the matrix elements are 

{n,q~q'\H^{x,y)\n,q) = f dk2V{q' ,k2)L^e-^^''"+'''^ei''''''e'''''' (30) 

Since the basis |n, q) form a complete set for the n-th level, it follows that the two operators 
Hf and Hf are equivalent in the one-band approximation if (28) and (29) are equivalent 
for any q and q' . Assume the equivalence of these two matrices, we can take inverse 
Fourier transforms of integrals with respect to q to get an identity between the integrands 
S(L^ + Li_i)exp(-(Q'2 + kj)/4)eMiQ'k2/2) and yL„exp(-(Q'2 + /e2)/4) exp(V/c2/2). 
After dividing out the common exponential part, it follows that Hi = Hf if and only if 

2L^V{kuk2) = {Ll + Ll_i)B{kuk2), ykuk2, (30), 

which is a continuous version of Eq.(27). The discussion in the previous subsection about 
the singularities in the Fourier components applies to this case as well. This relation might 
be helpful to the researchers working on problems of random magnetic field. 

VI. CONCLUSION 

We have calculated the energy spectra of 2D electrons in two dimensional periodic 
magnetic fields. The spectra are expressed in quasi-momenta that are good quantum 
numbers of the system. It is found that, the Landau levels are broadened in an oscillatory 
manner with respect to the band index. When the modulation of the magnetic field is 
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strong, dispersion curves are mixed but show avoided crossings. The magnetic lattice with 
triangular symmetry can be realized by using a superconductor film in a vortex state. A 
connection between the electric properties of the 2DEG and the magnetic properties of the 
vortex arrays is studied. It shows a possible way of getting information about the vortices 
through measurements of the bandwidth of the 2DEG. Finally, the connection between 
electric field modulation and magnetic field modulation is discussed. We find the effective 
electric potentials for weak 2D magnetic modulations. The effective potentials differ from 
band to band, and do not always exist. 
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Figure captions 

Fig. 1. (a), A plot of the magnetic field B{x,y) with triangular lattice symmerty in 
the 6-point approximation, (units of Bq) (b), The ground state energy of an electron in a 
homogeneous magnetic field modulated by the magnetic field in (a), (units of huJc) 

Fig. 2. Broadened Landau levels for a 2DEG in a homogeneous magnetic field mod- 
ulated by a square magnetic lattice. The perturbing field is B{x,y) = 0.5(cosx + cosy). 
Dashed lines are for the first order energy perturbation, and solid lines are for the exact 
energy spectrum. Only the lowest fifteen levels are shown here, (units of huc) 

Fig. 3. Exact energy spectrum for a 2DEG in a homogeneous magnetic field modulated 
by a triangular magnetic lattice in the 6-point approximation, 5(1, 0) = 1/3. 

Fig. 4. Exact energy spectrum for a 2DEG in a homogeneous magnetic field modulated 
by a triangular magnetic lattice in the 6-point approximation, 5(1,0) = 1. Many of the 
curves at the upper left and right are actually separated with higher resolution. 

Fig. 5. Hall conductivity at zero temperature. Due to the broadened Landau levels, 
the transition from one plateau to another is gradual. -B(l, 0) here is 0.2. 
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